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ABSTRACT 

Most of the massive elliptical galaxies in the universe stopped forming stars billions of years ago, even 
though plenty of hot gas remains available for star formation. Here we present compelling evidence 
indicating that quenching of star formation depends on both black-hole feedback and Type la super¬ 
nova heating. We analyze Chandra X-ray observations of ten massive ellipticals, five with extended, 
potentially star-forming multiphase gas and five single-phase ellipticals with no star formation. The 
ratio of cooling time to freefall time at 1-10 kpc in the multiphase galaxies is t coo \/tg ~ 10, indicating 
that precipitation-driven feedback limits cooling but does not eliminate condensation. In the same 
region of the single-phase galaxies, the radial profiles of gas entropy are consistent with a thermally 
stable (CooiAff > 20) supernova-driven outflow that sweeps stellar ejecta out of the galaxy. However, 
in one of those single-phase ellipticals (NGC 4261) we find t coo i/tg < 10 at < 300 pc. Notably, its jets 
are ~ 50 times more powerful than in the other nine ellipticals, in agreement with models indicating 
that precipitation near the black hole should switch its fueling mode from Bondi-like accretion to 
cold chaotic accretion. We conclude by hypothesizing that particularly strong black-hole outbursts 
can shut off star formation in massive elliptical galaxies by boosting the entropy of the hot gas and 
flipping the system into the supernova-sweeping state. 

Subject headings: galaxies: clusters: intracluster medium 


1. INTRODUCTION 


Abundant circumstantial evidence supports the hy¬ 
pothesis that black-hole feedback su ppresses star for¬ 
matio n in massive galaxies (e.g., McNamara & Nulsen 


2007 2012) but a deep mystery remains: How does ac¬ 
cretion fueling of the black-hole engine become closely 
coupled with the vast hot-gas atmosphere surrounding 
it? Precipitation-driven feedback models are providing 
particularly promising answers to this question. 

Numerical simulations show that cool clouds can pre¬ 
cipitate out of a galaxy’s hot-gas atmosphere via thermal 
instability if it is in a state of global thermal balance, 
with heating approximately equal to cooling (McCourt 


|et al. l 2012| ISharma et aT] 2012b| |Gaspar i et al. 2012[ see 


also Pizzolato & Soker (|2005[)). The critical criterion for 

precipitation depends on the ratio between the time t coo i 
required for gas at temperature T to radiate 3/cT/2 per 
particle and the free-fall time tg = ( 2r/g j 1 / 2 required for 
a dense cool cloud to fall from a radius r at the local grav¬ 
itational acceleration g. When the average f C ooiAff ratio 
is less than 10, then cooling is fast enough for some of 
the hot gas to condense into cold clouds and precipitate 
out of the hot medium. 

Precipitation itself plays an essential role in maintain¬ 
ing the required state of global thermal balance because 
it provides fuel for accretion. Recent numerical simula¬ 
tions of this feedback loop have been very encouraging, 
because they show that “chaotic cold accretion” of pre- 
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cipitating clouds can produce a black-hole fueling rate 
two orders of magnitude greater than the Bondi rate and 
can therefore generate a feedback response that brings 
the system into approximate balance at t rn oi / tg ~ 10 
flGaspari et al. 2013 201 4 Li & Bryan||2014a|b see also 
Pizzolato fe Soker! (2010)). 

We have recently presented measurements of galaxy- 
cluster atmospheres that strongly support this picture. 
The amount of multiphase gas in the cores of galaxy clus¬ 
ters steeply anticorrelates with min(t coo iAff), indicating 
that the precipitation rate, and therefore the feedback 
response, de pends very strongly on this parameter (Voit 
& Donahue 2015). Even more compelling is the finding 
that AooiAff ~ 10 is a global lower limit on this ratio at 
all radii and at all tempera tures amon g galaxy clusters 
in the A CCEPl ^ database ( Voit et aL|2014 |. 

Here we sho w, using Chand r a obse r vation s originally 
presented by | Werner et al (2012 2014), that the 

tcooiAff ~ 10 precipitation limit also applies to massive 
elliptical galaxies. Those observations revealed a close 
link between the presence of extended multiphase gas 
in ellipticals and the thermodynamic properties of their 
hot-gas atmospheres. This paper interprets that find¬ 
ing in terms of a model combining supernova-driven out¬ 
flows and precipitation-driven feedback. Section 2 briefly 
summarizes the key features of the Werner et al. ellip¬ 
ticals. Section 3 shows that Type la supernovae (SNIa) 
can drive outflows that sweep stellar ejecta out of ellip¬ 
ticals with velocity dispersion a v > 250kms _1 but can¬ 
not prevent precipitation within the central kiloparsec. 
Section 4 argues that outbursts of precipitation-driven 
black-hole feedback at < 1 kpc prevent t coo i from falling 
much below lOfff. Section 5 shows that the black hole’s 
jets are poorly coupled to the 1 10 kpc region but usually 

|2009[ ), 


ICavagnolo 


et 


al. 
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manage to maintain an isentropic core at < 0.5 kpc. Sec¬ 
tion 6 concludes by hypothesizing that black-hole feed¬ 
back shuts off star formation in massive ellipticals by 
switching on supernova sweeping. 


2. RADIAL PROFILES & PRECIPITATION 

Profiles of electron density (n e ), temperature (T), and 

entropy index (AT = fcTn e 2 ^ 3 ) for t en ma ssive elliptical 
galaxies from Werner et al. (2012 2014) are shown in 
Figure 1 as functions ot radius r = (1 kpc)r kpc . Dashed 
lines and blue symbols show the subset of five with mul¬ 
tiphase emission-line nebulae at r ~ 1 10 kpc. Solid 
lines and red or purple symbols show the subset of five 
without extended nebulae. The density panel shows 
that n e ~ (6 x 10~ 2 cm“ 3 )r kpc 3 is a good approxima¬ 
tion at 1 10 kpc for the single-phase ellipticals, whereas 
the multiphase ellipticals tend to have density profiles 
closer to n e oc r^ 1 . The temperature panel shows that 
the multiphase galaxies are systematically cooler than 
those without extended nebulae. In this sample, the tem¬ 
perature difference corresponds to a difference in veloc¬ 
ity dispersion: The five galaxies with extended nebulae 
have one-dimensional velocity dispersions in the range 
< t v = 217 -255kms _1 , and the five without extended neb¬ 
ulae have cr v = 263-336 km s -1 , according to the Hyper- 
Leda database^] 

Combining n e (r) and T(r) to obtain K{r) reveals that 
the hot atmospheres of multiphase ellipticals distinctly 
differ from those without extended nebulae. The single¬ 
phase elliptical galaxies follow the power law K(r) ~ 
(5 keVcm 2 )r kpc in the 1-10 kpc range, whereas the 
multiphase elliptical galaxies tend to follow K(r) ~ 

(3.5 keV cm 2 )r^ 3 . 

This latter power law corresponds to the pre cipitation 
threshold predicted by Sharma et al. ( 2012a) a n d con - 
firmed in central cluster galaxies by Voit et al. (2014). 
Setting t coo i = 10 tg « lOrcr” 1 yields 


-Kprecip(f) P3 (3.5 keV cm ) A 3 g_ 23 a 2 ^J r k p C , 

(1 ) 

where A = (3 x 10 23 erg cm 3 s 1 )A 3e _ 2 3 is the usual 
radiative cooling function for gas in collisional ioniza¬ 
tion equilibrium and depends on elemental abundances 
in this temperature range, with A 3e _ 2 3 ~ 1 for a solar- 
abundance plasma. Just as in central cluster galaxies, 
the entropy profiles of elliptical galaxies with extended 
multiphase gas at 1-10 kpc appear to be maintained 
by precipitation-driven feedback, strengthening the over¬ 
all case for this mechanism. However, something else 
must be happening at these radii in single-phase ellipti¬ 
cal galaxies. 


3. SUPERNOVA SWEEPING 

We propose that the difference between multiphase and 
single-phase ellipticals arises because SNIa are success¬ 
fully sweeping gas ejected by the old stellar population 
out of the single-phase ellipticals. Stars belonging to an 
old stellar population shed matter with a specific mass 
ejection rate 

a = a_i 9 x 10^ 19 s _1 ~ (3 x 10 -12 yr _1 ) a_ig , (2) 

2 http://leda.univ-lyonl.fr 


where a_ig is within a factor of 2 of unity but depends 
somewhat on the age of the stellar population and the 


low-mass tail of the stellar m ass function (e.g., Math- 
ews||1990 Conroy et al.||2014). The ejecta cannot linger 


because otherwise the galaxy would either contain more 
gas or would be forming stars at a rate e quivale nt t o the 
mass ejection rate (e.g., Mathews fe Brighenti|2003 Voit 


& Donahue 2011). Instead, the time-averaged outflow 
rate of stellar ejecta at radius r must be 


M(r) « (0.09 M 0 yr 3 ) a_i 9 er^so r kpc , (3) 

where we have assumed that the stellar mass-density pro¬ 
file is approximately oc r^ 2 with a one-dimensional veloc¬ 
ity dispersion a v = (250 km s _1 ) 02 . 50 , for which the stel¬ 
lar mass within radius r is ~ (1.4 x lO lo M 0 )cr| 5O r kpc . 
Given the power-law electron density profile shown in 
Figure 1, the average outflow velocity is then 

v « (Skms^ 1 ) 0-19 o-lso r k P 3 c , (4) 


implying that a steady outflow would be subsonic and 
therefore close to hydrostatic equilibrium. 

In order to illustrate how supernovae can produce 
an outflow with the power-law entropy profile observed 
in single-phase ellipticals, let us temporarily ignore ra¬ 
diative cooling. In that case the entropy equation for 
steadily outflowing gas can be expressed as 

din A' 2 H 5 ap * 

dt 3 nkT 3 pm p n ’ 


where H is the heating rate per unit volume, n is the 
number density of gas particles, pm p is the mean mass 
per gas particle, and p* is the local stellar mass den¬ 
sity. The first term accounts for the change in specific 
entropy due to heat input, and the second accounts for 
introduction of new particles into the system. 

Supernovae are not the only heat source. Mass ejected 
from old stars carries a specific kinetic energy per particle 
« 3p,TO p (7 2 /2 that becomes thermalized in the ambient 
medium as the ejecta merge with it (e.g., Conroy et al. 
2014). The specific energy ejected by the old stellar pop¬ 
ulation (3fcTjj/2) therefore separates into a supernova 
heating term (3/cTsn/2) and a kinetic energy term: 


= Hfmip « | [/cTsn + pm p al\ , (6) 

2 ap* 2 

Substituting for H in the entropy equation and multiply¬ 
ing by ru _1 then leads to 


din A' r T h 5 
d In r T 3 

where p* is the mean stellar mass density within r. In 
other words, the power-law slope of the radial entropy 
gradient depends primarily on the ratio of specific en¬ 
ergy per ejected particle to thermal energy per ambient 
particle. If that ratio exceeds 5/3, then the entropy gra¬ 
dient of the outflow will be positive, and if it does not, 
then the entropy gradient will be negative and convec- 
tively unstable. Furthermore, the density factor can be 
neglected for the stellar mass distribution we have as¬ 
sumed, in which p* « 3p*. 


3p* 


— 5 

P* 


(7) 
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Fig. 1. — Radial profiles of electron density n e (top panel), gas temperature T (middle panel), and entropy index K = kTn e 2//3 (bottom 
panel) in massive ellipticals with extended emission-line nebulae (dotted lines, blue symbols) and without extended nebulae (solid lines, 
red or purple symbols). A solid grey line in the top panel shows n e (r ) = (6 x 10 —2 cm -3 )r^ p 1 c 3 . A long-dashed grey line in the bottom 

panel shows K(r) = (5 keV cm 2 )rk pc . A solid pink line in the bottom panel shows K(r) = (3.5 keV cm 2 )r^, which corresponds to the 

precipitation threshold at £ C ooi/% ~ 10 for galaxies in this mass range. A short-dashed green line shows K{r) = (5keVcm 2 )r^, which 
corresponds to the unstable locus at which SNIa heating balances radiative cooling. 


This relationship conveniently allows us to infer the su¬ 
pernova heating rate from the slope of the entropy gradi¬ 
ent at 1-10 kpc. For the single-phase ellipticals, in which 
K oc r, we obtain 


o 

/cTsn ~ -kT — fjbm p (Ty ~ 2keV (8) 

o 

since they have kT ss 1.0 keV and a v ss 300kms _1 . In 
order to impart this much energy to the matter ejected 
by the old stellar population, the specific SNIa rate needs 
to be 


0.3 (100yr)" 1 (lO^Mo)" 1 ^ 1 ^ 


(9) 


where £sn = (10 51 erg)I? 5 i is the heat energy introduced 
per supernova into the ambient gas. 

This result is broadly consistent with direct ob serva¬ 
tions of S NIa rates in old stellar populations (e.g., Maoz 
et al.|20f21 and validates our outflow analysis. However, 
it is in tension with the iron abundances observed in the 


hot gas, which are approximately solar. An iron yield 
~ 0.7 Mg per supernova is equivalent to a mass fraction 
~ (7x lO -3 ).^ 1 relative to all the stellar ejecta, amount¬ 
ing to ~ 5-E.y times the solar abundance. Reconcilia¬ 
tion of the observed gas-phase iron abundances with the 
observed SNIa rate would therefore seem to require the 
newly-produced iron to be poorly mixed with the ambi¬ 
ent g as, at least within the central 10 kpc of the outflow 
(e.g., Tang et aL|2009 ). 

Now fet us return to the issue of radiative cooling. The 
cooling-free outflow solution for single-phase ellipticals 
must become invalid at small radii, because radiative 
cooling per unit volume scales as n 2 oc r~ 2 6 while heat¬ 
ing scales as p * oc r~ 2 . Defining n e ,SNbai = H 1 ^ 2 A -1 / 2 
and AgNbai = kTri e g(f b a| leads to the following expres¬ 
sion for the unstable locus of supernova-heating balance 
for a galaxy with kT « 0.8keV and <j v ss 300kms -1 : 


XsNbai ~ (5 keV cm 2 ) a_^ 3 A^ 3 23 r^ 3 . 


( 10 ) 
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This is shown with a green dashed line in the bottom 
panel of Figure 1. The power-law profiles of density 
and entropy observed in single-phase ellipticals do indeed 
break where they intersect this locus, indicating that su¬ 
pernova heating exceeds radiative cooling at larger radii, 
and that supernova sweeping alone cannot rid the inner 
~ 1 kpc of its stellar ejecta. 


4. BONDI ACCRETION & PRECIPITATION 


Inside of ~ 0.5 kpc, nine of the ten galaxies in our sam¬ 
ple become isentropic, with core entropy Kq ss 2 keV. 
These all have radio power ~ 10 38 ergs _1 . Adopt- 
ing the conversion from radio power to jet power from 


Cavagnolo et al. (2010) indicates a typical jet power 
~ 2 x 10"“ ergs -1 . The total power av ailable from Bondi 
accretion (see also Allen et al. 2006) can be estimated 
from the typical core entropy and black hole mass (Mbh) 
for these galaxies: 



Af Bo ndi c 2 ~ 4tt G 2 Mb H c 2 (/xto p ) 5/2 ( U ) 

~ (5 x 10 44 ergs" 1 ) JWfg H , 9 K^ /2 , ( 12 ) 


-3/2 


where Mbh ,9 = Mbh/10 9 M q and K 2 = A' 0 /2keVcm 2 . 
Dividing the jet power by this value gives an inferred 
jet efficiency ~ 10 —2-4 , meaning that the jets in these 
galaxies can plausibly be fueled by Bondi accretion onto 
the central black hole at the present time. 

The tenth galaxy is NGC 4261, shown with purple stars 
in Figure 1. It is the only one that drops below the pre¬ 
cipitation threshold within 0.5kpc (see Figure 2), and 
its radio luminosity (~ 3 x 10 4 O ergs _1 ) corresponds to 
a jet power ~ 10 44 ergs _1 , almost two orders of magni¬ 
tude greater than in the other galaxies. From the cen tral 
black-hole mass of this galaxy (5 x 10 8 M©, Kormendy 
fc Ho||2013| ) and its entropy in the innermost bin, we in- 
fer a total Bondi power ~ 3.5 x 10 44 erg s _1 , which would 
imply an implausibly large jet efficiency ~ 0.3. We there¬ 
fore suggest that NGC 4261 may be a case in which cold 
chaotic accretion triggered by precipitation of cold clouds 
has temporarily boosted the accretion rate far above th e 
Bondi level, as predicted by Gaspari et al. (2013 2014). 

In all cases, the jet power greatly exceeds the 
10 41 ergs _1 X-ray luminosity from < 1 kpc, implying 
that the jets are poorly coupled to the isentropic core and 
deposit most of their energy at larger radii. This find¬ 
ing agrees with the radio morphologies of the jets, which 
extend well beyond the central kpc. Tapping a small 
percentage of their power would compensate for cooling 
within ~ 1 kpc, and the flatness of the inner entropy pro¬ 
files suggests that some combination of turbulent mixing 
and convection distributes that heat energy. Further¬ 
more, if the jet ever fails to couple to the isentropic core, 
the core will fill up with stellar ejecta, thereby lower¬ 
ing Kq until the system reaches the precipitation limit. 
According to models of cold chaotic accretion, the jet 
power should then rise by two orders of magnitude, as ob¬ 
served in NGC 4261. In this sense, precipitation-driven 
feedback is like a backup system which ensures that the 
central cooling time remains > 10 tff, should either super¬ 
nova sweeping or Bondi accretion ever fail to keep stellar 
ejecta from accumulating within the galaxy. 


Fig. 2. — The precipitation criterion t coo i/tff as a function of 
radius in single-phase ellipticals (top panel) and multiphase ellip¬ 
ticals (bottom panel). Single-phase ellipticals all remain above the 
(pink) precipitation zone at 5 < t coo \/tg < 20, with one exception: 
NGC 4261 (purple stars) drops below the precipitation threshold 
at. Aooi/hi ~ 10 (magenta line) within the central kpc. In contrast, 
the multiphase ellipticals track the precipitation threshold in the 
1—10 kpc region where multiphase gas is found. Notably, the jets 
from NGC 4261 are two orders of magnitude more powerful than 
those from the other nine galaxies in this sample, indicating that 
chaotic cold accretion of precipitating clouds, not Bondi accretion, 
is fueling the black-hole engine in this system. 


5. JET COUPLING 


A simple scaling argument suggests that the jets couple 
poorly to the inner regions of ellipticals because of their 
power. In an isothermal potential, gas at the precipi¬ 
tation threshold has K oc r 2//3 and \dlnT/dlnr\ <C 1, 


implying n e oc r 

1 and kT ss 2/tm p cx 2 « 0.8 cr | 50 

keV. 

Setting t c ooi ~ 10tff therefore gives 


K(r) » 

i ( 3 ' 6keVcm2 ) A 3e-23 r kpc 

(13) 

n e (r) « 

i (0.1 cm 3 ) A 3e-23 *250 r kpc 

(14) 

Lx(< r) r 

* (10 41 erg s" 1 ) A ^ e 1 _ 23 crf 50 r kpc 

(15) 


where Lx is the X-ray luminosity. Following |Voit fe Don-| 


ahue (2005|, we estimate the shock velocity u s h driven 
through gas with this density profile by a jet of power 
Pjet = ( 10 42 ergs _ 1 )P 4 2 by setting m p n e r 3 v^ h ~ rv~^P iet 
and obtaining u s h ~ ( 800 kms- 1 )P 42 ^ 3 A 3 ^ 23 cr^ 5 Qrj( p 1 / 3 . 
For this sh ock ve locity, the entropy jump condition from 


Voit et al. (2003) yields an entropy increment 


AK „ H m P V sh 

~ 3(4n e ) 2 / 3 

~ (2- 4 keV cm2 ) P 4 2 2 /3 A 3e—23 *250 


(16) 


that is independent of radius. 

In other words, jet outbursts of the power observed in 
nine of these ten systems suffice only to boost Kq to the 
observed level and have little effect on the ~ 1-10 kpc 
region. That is because they deposit the bulk of their 
power at greater radii. Most of it is probably thermal- 
ized beyond the radius r t h ~ (10 kpc)P42cr^5gA3 e _23 at 
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which jet-driven shocks become subsonic. This is also 
the radius at which the typical jet power is comparable 
to the X-ray luminosity in the precipitating configura¬ 
tion, which is probably not a coincidence. Furthermore, 
the excellent agreement of single-phase ellipticals with 
the power-law supernova sweeping model, along with the 
lack of scatter in their K(r) profiles at < 10 kpc, sup¬ 
ports the idea that most of their jet power is deposited 
farther out. Even NGC 4261 agrees well with the oth¬ 
ers at 1-10 kpc, implying that the vast majority of its 
~ 10 44 ergs -1 jet power passes straight through the 1-10 
kpc region with very little dissipation. 


6. THE BLACK-HOLE FEEDBACK VALVE 

We conclude with a hypothesis for the cessation of star 
formation known as quenching. It would appear from the 
data presented in §2 that black-hole feedback alone can¬ 
not completely quench massive elliptical galaxies because 
some have persistent multiphase gas. Instead, massive 
elliptical galaxies separate into two groups. In single¬ 
phase ellipticals, SNIa keep stellar ejecta from accumu¬ 
lating by sweeping it out of the central 1 10 kpc (§3). 
In multiphase ellipticals, supernova heating cannot quite 
overcome radiative cooling in the 1 TO kpc region. These 
galaxies therefore fall into a precipitating state in which 
black-hole feedback maintains f coo i ~ lOtff and can can¬ 
not entirely shut off star formation (§4). Jet power can 
fluctuate between ~ 10 42 ergs -1 (Bondi-like accretion) 
and ~ 10 44 ergs -1 (cold chaotic accretion) without dis¬ 
rupting the 1 10 kpc region because the jets deposit most 
of their energy beyond ~ 10 kpc (§5). 

Why then are the five ellipticals with o v > 260 km s -1 
fully quenched while low-level star formation persists in 
the five multiphase ellipticals with lower velocity disper¬ 
sions? Some insight can be gained by comparing A' prec j p 
with i\ SNbal . In galaxies with A precip > A' SNbal , the 
interstellar medium can be thermally unstable even if 
supernova heating exceeds radiative cooling. We there¬ 
fore expect those galaxies to experience multiphase cir¬ 


culation in which SNIa drive hot gas out of the cen¬ 
tral regions where the stars are but cool clouds pre¬ 
cipitate back in from the margins. In contrast, galax¬ 
ies with AgNbai > A^rip can drive outflows capa¬ 
ble of shutting off precipitation, by raising f coo i ss 

(5Gyr)T“y /2 A 3e _23(A7l00keVcm 2 ) 3 / 2 high enough to 
prevent the outflowing gas from cooling after it leaves 
the galaxy. 

The condition for avoiding potentially star-forming 
multiphase gas is then 


K 

1 v precip 

A'SNbal 


0.6 a 


1/3 

-19 


A 


1/3 
3e—23 


—4/3 

7 250 


<1 , 
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implying that galaxies with a v < 

^ 3 e -23 should be embedded within a precipitating cir- 
cumgalactic medium. Intriguingly, this result implies 
that quenching depends more directly on the concentra¬ 
tion of stellar mass toward the center of a galaxy than 
on the total stellar ma ss itself, in qualitative agreement 
with observations (e.g., |Fang et ah [2013). It also implies 
that quenching requires a greater value of a v at earlier 
times, when the stellar population is younger and a is 
greater. However, the bottom panel of Figure 1 implies 
that something more is needed for quenching to be com¬ 
plete. 

We hypothesize that black-hole feedback outbursts are 
also required to flip massive ellipticals with A'gNbai > 
Aprecip from the precipitation locus to a supernova¬ 
sweeping state. This “black-hole feedback valve” can 
accomplish the task with an outburst powerful enough 
to produce an entropy increase AK > AsNbai — Ap re cip. 
It will be interesting to test this quenching mechanism 
and see how it plays out in simulated elliptical galaxies. 
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